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Abstract 

Many studies of anomalous microwave emission (AME) have computed an AME 
emissivity to compare the strength of the AME detected in different regions. Such a 
value is usually defined as the ratio between the intensity of the AME at 1 cm and the 
thermal dust emission at 100 fim. However, as studies of Galactic dust emission have 
shown, the intensity of the thermal dust emission at 100 fim is strongly dependent on 
the dust temperature, which has severe implications for the AME emissivity defined in 
this way. In this work, we illustrate and quantify this effect and find that the AME 
emissivity decreases by a factor of 11.1 between dust temperatures of 20 and 30 K. 
We, therefore, conclude that computing the AME emissivity relative to the 100 //m 
emission does not allow for accurate comparisons between the AME observed in different 
environments. With this in mind, we investigate the use of other tracers of the dust 
emission with which to compute the AME emissivity and we ultimately conclude that, 
despite the difficulty in deriving its value, the column density of the dust would be the 
most suitable quantity with which to compute the AME emissivity. 
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1 Introduction 



In recent years there has been growing evidence for the existence of a new component of 
microwave emission present in the interstellar medium (ISM). This emission appears to be 
spatially correlated with the dust in the ISM, although it produces an excess of emission 
with respect to the predicted thermal, vibrational dust emission at these frequencies (e.g. 



Leitch et al. , 1997), and as such this component is often described as anomalous. There have 



been only a handful of detections of this anomalous microwave emission (AME) component 



originating from both the diffuse ISM at mid-to-high latitudes (e.g. Banday et al. , 2003 



Collaboration, 2011b 



Scaife et al. 



Tibbs et al. 



Davies et al. 2006, Ghosh et al. , 2012) and specific Galactic sources such as the Perseus and 



Tibbs et al. 



2009), and a variety o 



p Oph molecular clouds (Watson et al. , 2005, Casassus et al. 2008, Tibbs et al. 2010; Planck 



2012b), the dark clouds LDN1622 (Finkbeiner et al., 2002 



Casassus et al. 2006), LDN1621 (Dickinson et al. 2010) and LDN1111 (AMI Consortium: 



HII regions (Dickinson et al. , 2007, Todorovic et al. , 2010 



2012a). Additionally, the first detection of AME in a star forming region in the 



external galaxy, NGC6946, was reported (Murphy et al. 2010, AMI Consortium: Scaife et 
al. 2010), confirming that this emission mechanism is truly ubiquitous. 

Observations have shown that the AME occurs in the frequency range ~ 10 - 100 GHz 
and is highly correlated with the mid-infrared (IR) dust emission. It is this correlation 



with the mid-IR emission that led Draine & Lazarian (1998) to postulate their model of 



spinning dust emission. In this model, the observed excess emission is produced by the very 
smallest dust grains (Very Small Grains; VSGs or Polycyclic Aromatic Hydrocarbons; PAHs), 
characterised by an electric dipole moment, rotating rapidly, resulting in the production of 
electric dipole radiation. The model predicts a well defined peaked spectrum, rising below 
~ 20 - 30 GHz and then falling off at higher frequencies as expected from a Boltzmann cutoff 



in grain rotation frequencies. The Draine & Lazarian (1998) spinning dust model has been 



Verstraete, 2010, Hoang et al. , 2011; Silsbee et al. 2011), and current spinning dust models 



2009, 


Hoang et al. 


2010; 


Ysard & 



incorporate a variety of grain rotational excitation and damping processes: collisions with 
neutral and ionised gas particles, plasma drag (the interaction between the electric field of 
ions and the electric dipole moment of the dust grains), absorption and emission of a photon, 
the photoelectric effect, microwave emission and the formation of H2 on the grain surface. 
Some of the more recent models incorporate additional features including dust grains that 



are not only spinning about their axis of greatest inertia (Hoang et al. , 2010; Silsbee et al. 



2011) and dust grains of irregular shape (Hoang et al. 2011) 



Understanding how the AME varies between different phases of the ISM is extremely 
important in improving our understanding of the AME. To help compare the strength of 
AME detected in one region to another, some authors have computed an AME emissivity. 
Given the strong association between the AME and the dust grains in the ISM, the AME 
emissivity is usually defined as the ratio between the antenna temperature of the AME at 
wavelengths of 1 cm and the surface brightness of the thermal dust emission at 100 /zm. 
This quantity, first calculated as a cross-correlation coefficient over large areas of sky (e.g. 
Kogut et al. 1996 de Oliveira-Costa et al. , 1997), has now been computed for a variety 



of individual regions in which AME has been detected (see Table Q. However, the AME 
emissivity defined in this manner is a highly biased method for making these comparisons 
because the surface brightness at 100 fim is significantly dependent on the dust temperature. 
In this work we aim to (1) illustrate the weakness of such a definition and (2) present 
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and discuss possible alternatives. In Section [2] we point out why the commonly used AME 
emissivity definition is flawed and estimate the amplitude of the bias it introduces when 
comparing regions with different dust temperatures. In Section [3] we discuss more suitable 
definitions of the AME emissivity including the use of gas and dust column densities. Finally, 
we present our conclusions in Section [4} 



Table 1: Sample of AME emissivity values for AME detections in a variety of different 
environments computed relative to the intensity of the 100 /xm emission. 



Source 



AME Emissivity 
OK (MJy sr- 1 )- 1 ) 



Reference 



Hll Regions 

6 southern Hi I regions 

9 northern Hi I regions 

Pleiades 

RCW49 

RCW175 

High Latitudes 

15 regions WMAP 
All-sky WMAP 



3.3 ± 1.7 
3.9 ± 0.8 
2.01 ± 0.09 
13.6 ± 4.2 
14.2 ± 2.7 



11.2 ± 1.5 
10.9 ± 1.1 



Dickinson et al. 


(2007 


') 


Todorovic et al. 


(201C 


)) 


Genova-Santos et al. 


(2011) 


Dickinson e 


al. 


(2007 


') 


Tibbs et al. 


(2012a 


) 



Davies et al. 



Davies et al. 



(2006) 



(2006) 



Molecular Clouds 

Perseus 
Perseus Al 
Perseus A2 
Perseus A3 
Perseus B 
Perseus C 



15.7 ± 0.3 
2.8 ± 0.7 
16.4 ± 4.1 

12.8 ± 6.1 
13.2 ± 3.6 
13.0 ± 3.2 



Watson et al. 



Tibbs et al. 



Tibbs et al. 



Tibbs et al. 



Tibbs et al. 



Tibbs et al. 



(2010) 



(2010) 



(2005) 



(2010) 



(2010) 



(2010) 



Dark Clouds 

LDN1621 
LDN1622 



18.1 ± 4.4 
21.3 ± 0.6 



Dickinson et a 



C 



asassus 



et al. 



(2010) 



2006) 



2 The AME Emissivity Defined Relative to the 100 fim 
Emission 

The concept of an AME emissivity is particularly useful as it provides a normalisation of the 
strength of the AME detected in a variety of regions, and hence it facilitates a comparison 
between such AME detections. In Table [T] we list AME emissivities defined relative to the 
100 jum emission from the literature for a variety of different phases of the ISM, including 
Hll regions, molecular clouds, dark clouds and diffuse emission at intermediate latitudes. 
Looking at the values in Table [T] it is possible to identify the observed range of values of the 
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AME emissivity. A typical value for the AME emissivity is of the order 10 /zK (MJy sr x ) 1 
with a range from ~ 3 to 25 fiK (MJy sr -1 ) -1 (Banday et al. , 2003, Davies et al. 2006). 



In terms of flux density, this corresponds to approximately 1 Jy of AME at 1 cm for every 
3000 Jy at 100 /im. 



Observations at far-IR and submillimeter wavelengths (e.g. Reach et al. 1995 Boulanger 



et al. , 1996) have revealed that the spectrum of the thermal dust emission is well represented 



by a modified black body function of the form 



S u oc 



2hv 3+f> 



1 



(1) 



c 2 exp(hv/kT dust ) - 1 

where (3 is the dust emissivity spectral index and Td us t is the equilibrium temperature of 
the big dust grains. The thermal dust emission observed at such far-IR and submillimeter 
wavelengths is produced by big dust grains absorbing ultra-violet photons from the exciting 
radiation field of the environment and re-radiating this emission thermally. It is this balance 
between absorption and emission that produces the well-defined spectrum of thermal dust 
emission (see Figure [T| and keeps the big grains in thermal equilibrium. In the ISM, the 



dust temperature is known to vary from one environment to an other (Planck Collaboration 



2011a) with dust temperatures ranging from ~ 15 K in dense molecular clouds to 25 - 30 K in 



Hll regions. Additionally, the dust temperature has been observed to vary on small angular 
scales and within regions in which AME has been detected (e.g. between 15 and 25 K in the 
Perseus molecular cloud; Tibbs et al. , 2011 ). Since (5 is dependent on the physical properties 



of the dust grains, it is also expected to vary between different environments: Pupae et 



al. (2003) found f3 to vary from 0.8 to 2.4 in a range of Galactic environments with dust 



temperatures between 11 and 80 K. Therefore, given that the emission at 100 /im is an 
observed quantity that depends strongly Td ust and /3, and since both Td us t and (3 are known 
to vary in the different phases of the ISM, the 100 /im emission is not a suitable quantity 
with which to estimate the AME emissivity. 

To quantify the effect of dust temperature variations on the AME emissivity, we consider 
a dust emissivity spectral index of 1.8, which is the value derived for the Galactic solar 



neighbourhood (Planck Collaboration, 2011a), and two dust temperatures (20 and 30 K) 



We then calculate that for a given intensity at 1 cm, the AME emissivity at 20 K is 11.1 
times higher than at 30 K. This effect is substantial and could by itself, without any intrinsic 
variation of the AME intensity, explain the variations of the AME emissivities listed in 
Table [TJ This clearly demonstrates the bias in the AME emissivity when computed relative 
to the 100 fim emission. With this in mind, we can now revisit the AME emissivity values 
listed in Table [TJ For example, we focus on the five regions of AME detected in the Perseus 
molecular cloud (Al, A2, A3, B and C) by Tibbs et al. (2010). It is apparent that region 
Al is much less emissive (by a factor of ~ 4 - 5) than the other four regions. However, 
as discussed by Tibbs et al. (2010, 2011), the physical environment of region Al, which 
corresponds to the open cluster IC348, is very different from the other four regions. In fact, 
simply because of the difference in dust temperature between IC348 (~ 22 K) and the other 
regions (~ 18 K), we estimate that the AME emissivity is reduced by a factor of 4.3 in region 
Al relative to the other regions. Therefore, the inconsistency between region Al and the 
other regions can be solely accounted for by dust temperature variations, illustrating the 
impact of the dust temperature in biasing the AME emissivity computed using the 100 /im 
emission. 
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Figure 1: The modified black body spectrum representing the thermal dust emission for a 
range of dust temperatures from 15 to 45 K with a fixed dust emissivity spectral index of 1.8. 
This illustrates how the 100 /im emission (vertical dashed line) increases with increasing dust 
temperature, and hence how the AME emissivity defined relative to the 100 /im emission is 
biased by changes in dust temperature. 



3 Redefining the AME Emissivity 



In Section [2] we have shown that the commonly used definition of the AME emissivity 
computed relative to the 100 /im emission is flawed. Therefore, an unbiased definition of the 
AME emissivity is required to help quantify the AME variations. 



As part of their investigation of the diffuse emission at intermediate latitudes, Davies et 



al. (2006) computed the AME emissivity relative to the dust emission at 94 GHz rather than 



100 /mi. Since 94 GHz is further from the peak of the thermal dust emission (see Figure [T 
the emission at 94 GHz is much less dependent on the dust temperature than the emission 
at 100 /im. However, even using the intensity of the 94 GHz emission, there is still an effect 
due to the dust temperature, as the AME emissivity decreases by a factor of 1.6 when the 
dust temperature increases from 20 to 30 K. 

Rather than using the intensity of the dust emission to compute an AME emissivity, other 
authors (e.g. Finkbeiner et al. 2002; Vidal et al. 2011) have computed an AME emissivity 



relative to the hydrogen column density, Nr. This definition has several advantages. First, it 
results in an AME emissivity that is more directly suited to comparing with the theoretical 
values as the current spinning dust theories compute an emissivity in units of Jy sr -1 cm 2 
per H atom. Using Nh also mitigates the effects of dust temperature, hence allowing a much 
more impartial comparison between measurements of AME in different environments. To 



illustrate this, we use the Nh map of the Perseus molecular cloud from Tibbs et al. (2011) 



which is based on a visual extinction map, and compute the AME emissivity for the five 
regions of AME in the Perseus molecular cloud (see Table [2j. When computing the AME 
emissivity relative to Nh, we find that region Al is now consistent (within 2 - 3a) with 
the four other regions in the Perseus molecular cloud, which leads to a completely different 
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interpretation. We now want to explore how this new definition can help us to perform 
a less biased comparison between AME detections in different environments. |Vidal et al. 



(2011) compared the AME emissivity of six regions covering a range of column densities 
from Galactic cirrus to dense dark clouds (see Table [2]). Adding these six measurements to 
the five we computed for the Perseus molecular cloud, gives a sample of 11 AME emissivities 
based on Nh- It is apparent that the AME emissivities for the Perseus molecular cloud are 
consistent with the values computed for other AME detections. It is also noticeable that 
there is less scatter in the AME emissivities listed in Table [2] compared to those listed in 
Table HJ 



Table 2: AME emissivity values for a variety of phases of the ISM computed relative to the 
hydrogen column density. 



Source 


N H 


AME Emissivity 


Reference 




(10 22 H cm" 2 ) 


(10" 18 Jy sr- 1 cm 2 /H) 




Perseus Al 


0.82 ± 0.19 


2.2 ± 0.5 


This work 


Perseus A2 


1.19 ± 0.27 


2.7 ± 0.6 


This work 


Perseus A3 


1.28 ± 0.29 


1.2 ± 0.3 


This work 


Perseus B 


0.79 ± 0.18 


3.4 ± 0.8 


This work 


Perseus C 


0.74 ± 0.17 


3.6 ± 0.8 


This work 



Cirrus 
C Oph 
LDN1780 
LDN1622 
p Oph 
M78 



0.15 ± 0.07 
0.22 ± 0.02 
0.45 ± 0.04 
1.50 ± 0.15 
5.00 ± 0.50 
22.80 ± 0.23 



4.6 ± 2.0 

4.1 ± 0.6 
3.5 ± 0.4 
2.0 ± 0.2 

3.2 ± 0.5 
0.9 ± 0.1 



Leitch et al. 



Vidal et al. 



Vidal et al. 



Casassus et 



(1997) 



2011) 



(2011) 



al. 



C 



asassus 



et al. 



(2006) 



Castellanos et a 



(2008) 



(2011) 



Following the example of Vidal et al. (2011), we plotted the AME emissivities listed 
in Table [2] as a function of Nh- This plot is displayed in Figure |2j We fitted the data 
with a power-law and found that the best fitting spectral index is —0.31 ± 0.03. Although 
there is some scatter in the data, the results hint that the AME emissivity is decreasing 
with increasing Nh, suggesting that the AME becomes less emissive as the column density 
increases. This is consistent with the AME being due to the smallest dust grains, as it is 
known that as Nh increases, the abundance of the smallest dust grains decreases due to dust 
grain coagulation (e.g. Stepnik et al. 2003, Flagey et al. , 2009). The idea of the AME being 
associated with the smallest dust grains is additional support for the AME being due to 
spinning dust. 

Although we have shown that defining the AME emissivity using Nh removes the bias in- 
troduced by the dust temperature and allows for a much more impartial comparison of AME 
between different environments, estimating Nh, including molecular, atomic and ionised 
forms, is technically challenging. There are three methods used to estimate Nh: near-IR 
extinction mapping; far-IR thermal dust emission; and the use of gas tracers. However, 
each method has its limitations. For example, near-IR extinction is only reliable if there 
are a sufficient number of background sources available, thermal dust emission estimates 
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Figure 2: AME emissivity computed using N H for the AME detections listed in Table [2] 
plotted as a function of Nr. Also plotted is the fit to the data (dashed line), which we 
estimate to have a spectral index of —0.31 ± 0.03. 



are affected by variations in the dust opacity, while observations of gas tracers, although 
they do provide kinematical information, are limited by the critical density of the observed 
molecule. A detailed analysis performed by Goodman et al. (2009) clearly illustrates the 
various uncertainties and limitations in estimating Nh using the three separate methods. 

In addition to the difficulties in measuring Nh, the AME is believed to be associated with 
interstellar dust grains, and hence it would make more sense to normalise the strength of the 
AME with a physical property of the dust. For example, the column density of dust, Nd us t, 
would be a useful quantity as it has the same advantages as using Nh with the addition 
that it is a direct tracer of the dust. However, Nd us t is even more difficult to estimate than 
N H . It is possible to assume a given M dust /M H and convert from N H to N dust , however 
this assumes that we know Mdust/Mn- Globally, within our own Galaxy, it is known that 
M dust /M H is 



0.01 



e.g. 



Draine & Li, 2007), but AME has recently been observed in external 



galaxies (e.g. Murphy et al. 2010 Planck Collaboration 2011c) where it is known that the 
metallicity, and hence Mdust/Mn, differs from our own Galaxy (e.g. Draine et al. , 2007). To 
mitigate the affects introduced by variations in Mdust/Mn, Ndust may be computed by fitting 
the far-IR SED with a dust model which incorporates the physical properties of the dust 
grains. However, to date, no such modelling has yet been performed systematically for the 
environments in which AME has been detected. 



4 Conclusions 

Having discussed various definitions of the AME emissivity, it is clear that using the 100 /im 
emission introduces a bias due to the effect of dust temperature, and that using Nh allows 
for more accurate comparisons. However, given the association between the AME and the 
interstellar dust, an intrinsic property of the dust such as Nd us t ; would represent the best 
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quantity with which to compute the AME emissivity. Therefore, until such a time as a new 
AME emissivity is defined, we stress that care should be taken before using it to directly 
compare the strength of AME in different environments. 
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